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INTRODUCTION 


The acoustic admittance of practical liners may be a function of the aeroacous- 
tic environment in which they are located. This fact necessitates that their proper- 
ties be evaluated in such an environment. This is accomplished in the laboratory 
with a grazing-flow impedance tube (or rectangular duct). In such a tube, the acous- 
tic material is aligned so that the sound and mean flow graze over the surface of the 
material in a controlled manner to simulate the environment typical of that found in 
a real jet engine. Since admittance is defined as the ratio of normal particle 
velocity to pressure at the wall surface, the admittance could be determined experi- 
mentally by direct measurements of normal particle velocity and pressure. In prac- 
tice, measurements of normal velocity in the presence of grazing flow are beyond 
present measurement techniques. Instead, an indirect approach must be used which 
depends only on the measurement of acoustic pressure at selected locations in the 
test configuration. These measurements are then used as input to an analytical pro- 
gram which determines the admittance based upon these measured data. 

To date, two methods have been presented for determining the acoustic admittance 
in grazing flow from measured values of the axial propagation constant and cross 
modes. The first method is that described by Mungur and Gladwell (ref. 1). This 
method is applicable to general flows and is based upon a Runge-Kutta integration 
technique. The second method is described by Armstrong et al. (ref. 2). Armstrong's 
method is restricted to thin boundary-layer flows. This method uses an asymptotic 
expansion to derive a boundary condition that is applicable at the boundaries of the 
uniform region, in which an exact solution is possible. 

Both methods presented in references 1 and 2 are based upon the solution to an 
ordinary differential equation and are restricted to mean flows with transverse gra- 
dients in one direction only. However, grazing-flow duct facilities commonly employ 
ducts with a relatively small rectangular cross section in which the grazing flow 
possesses gradients in both transverse directions of the impedance tube. Such is the 
case here in the flow impedance test laboratory at the Langley Research Center. 

Thus, the present effort was motivated by the need to account for the more realistic 
flow environment in laboratory flow-impedance tubes such as those in the Langley 
facility. 

This paper describes a new method for determining the admittance of an acoustic 
material in a grazing-flow impedance tube (or duct). The method is developed explic- 
itly for rectangular ducts, although it is applicable to other duct geometries as 
well. The approach is to obtain an estimate of the axial propagation constant, 
cross-mode order, and mean-flow profile by means of acoustic and aerodynamic measure- 
ments. These experimentally determined characteristics of the aeroacoustic field 
are then input into an analytical program from which the admittance is calculated. 
Unlike the methods employed in references 1 and 2, the approach presented in this 
research is extendable to include mean flows with gradients in both transverse direc- 
tions of an impedance tube. This paper gives a detailed description of the method 
employed when the mean flow is restricted to gradients in only one transverse direc- 
tion. A formulation of the method to include mean-flow gradients in both transverse 
directions will be reserved for a future paper. 



The present paper is divided into three sections. The first section presents 
the governing differential equation, boundary conditions, and a description of the 
measured data necessary for a determination of the admittance. The second section 
describes the numerical method employed to determine the admittance. The paper 
closes with a results section (third section) in which admittance values determined 
by this method are compared both with exact values that can be obtained for a con- 
stant mean-flow profile and with results from a Runge-Kutta integration technique for 
cases involving a one-dimensional boundary layer. Such comparisons serve as confi- 
dence builders giving credibility to the method. 


SYMBOLS 


[A(p M JMBMdH 

\ global matrices 
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[B] , [5] 


[A q ] 


stiffness matrix for element q 
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c 

E(y) 

F(y) 
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h 

i 

K 

K 

x 

A# 

K x 

L 

M(y ) 


dimension of line element 
ambient speed of sound 
residual error 

one-dimensional acoustic pressure eigenfunction 
respective values of functions F(y) and M(y) at node I 
element shape function 
duct height 

free-space wave number 

dimensional axial propagation constant 
dimensionless axial propagation constant 
duct width 

Mach number of mean-flow profile 
centerline Mach number 
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n 

P(y,z) 

t 
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P w,u' P w,A 
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Subscripts: 


cross -mode order 

total number of locations in transverse direction of impedance tube 
boundary-condition parameter 

acoustic pressure eigenfunction in impedance -tube cross section 
time 

rectangular coordinates 

local coordinate for element 

value of global coordinate y at node I 

specific acoustic admittance of upper and lower walls, respectively, 
impedance tube 

eigenvalue 

global vector of unknowns 

vector of unknown nodal values for element q 
angular frequency 


I,J 

N 


integers 

number 


Superscripts: 

n integer (1 or 2) 

q element number 

Tr transpose 

Mathematical notation: 

[ ] denotes matrix 

{ } denotes vector 

• first derivative 


M 


second derivative 



PROBLEM FORMULATION 


Figure 1 depicts the cross-section geometry of the grazing-flow impedance-tube 
(rectangular-duct) test section and its coordinate system. The mathematical 



Figure 1.- Grazing-flow impedance -tube test section and coordinate system. 


development assumes that the mean flow in the test section is axial and fully devel- 
oped so that it is independent of the axial coordinate x. The two sidewalls of the 
tube are assumed rigid, whereas the upper and lower walls are lined with acoustic 
material. The specific acoustic admittances of the upper and lower walls are P w 
and p w respectively. All four walls in the figure are assumed to be of infinite 
extent in the axial direction (infinite-length assumption), and sound waves are 
assumed to travel in the axial direction. Further, the mean flow is assumed to have 
an axial component only and flow gradients are assumed to exist only in the 
y-di recti on of the tube. 

Upon taking into account the fundamental equations of fluid dynamics (i.e., the 
momentum equations and the continuity equation) and assuming that the acoustic pres- 

-i(K x-cot) 

sure field is of the form P(y,z) e X , the following elliptic partial 

differential equation for the pressure eigenfunction P(y,z) is obtained (ref. 3): 


P + P 

yy zz 


+ 


2K M P 
x y y 


1 - MK 




( 1 ) 


where 


K = v* 
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In equation (1) the function M(y) is the Mach number of the mean-flow profile, 
is the dimensionless axial propagation constant in the tube, and K is the free- 
space wave number u)/c. 

Boundary conditions associated with equation (1) are expressed in the form 


P z (y,o) = o 

(3) 


(4) 


(5) 


Continuity of acoustic-particle displacement at the lined walls is obtained by set- 
ting the integer n to 2, whereas continuity of acoustic-particle velocity is 
obtained by setting n to 1 . In recent years, continuity of acoustic-particle 
displacement has been more universally accepted as the correct boundary condition; 
however, realistic mean-flow profiles satisfy the condition of no slip at the bound- 
aries and render a zero mean flow there. Thus, for realistic mean flows, both forms 
of the boundary condition coalesce to the same expression at each wall. However, to 
preserve the generality, both forms of the boundary condition will be incorporated in 
this paper by carrying n throughout the analysis as a parameter. 

The solution to equations (1) through (5) can be expressed in the form 


P z (y,L) = 0 


p y (0,z) = 


1 - M( 0 ) K 


~ n 

P(0, 

x 


z) 


P y (h,z) = 


1 - M(h) K P(h,z) 


H n 
x 


P(y,z) = F(y) cos(mitz/L) (m = 0, 1, 2, 3, ...) 


( 6 ) 


where the eigenfunction F(y) satisfies the ordinary differential equation 


F" 


+ 


2K M'F' 
x 






+ (irm/KL) 
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>F 


= 0 


(7) 


and the boundary conditions 


F' (0) 
F' (h) 


-iKp 


w,A 


1 - M(0) K 


n 


F( 0 ) 


iK{3 


w,u 


1 - M(h) K 


n 


F(h) 


( 8 ) 

(9) 


5 



Equations (6) through (9) are easily derived by finding the general solution to equa- 
tion (1) by separation of variables and then applying the boundary conditions (equa- 
tions (2) through (5)). 

Equations (7) through (9) constitute a boundary-value problem for the acoustic 
pressure eigenfunction F(y). The wave number K, mean-flow profile M(y), and 
boundary-condition parameter n in these equations are specified and are treated as 
known# ^The four remaining parameters of this problem are the axial propagation con- 
stant k , the cross-mode order m, and the admittances of the wall lining p w u 
and P w Only three of these parameters are independent and can be specified. 

The fourth parameter is then determined so that the governing differential equation 
and boundary conditions are satisfied. Generally, when qiven the admittances 
and p w ^ and cross-mode order m, a corresponding set of axial propagation con- f 

stants K x are determined. This paper focuses on the inverse problem. Thus, the 
admittance P w ^ will be specified, K x and m are assumed to be obtained from 
acoustic measurements, and the admittance 8 t _ will be determined. Although an 

W f XJi 

r* 

infinite set of axial propagation constants K x exist for each cross-mode order m, 
it is assumed that one member of this set can be measured. 

Upon obtaining the value of the necessary parameters, the boundary- value problem 
posed by equations (7) through (9) can be solved. However, equation (7) is an ordi- 
nary differential equation with variable coefficients, and its solution can be put in 
terms of known functions only for some special cases of M. Therefore, these cases 
are not useful for general application. Thus, the admittance must be determined 
numerically for mean flows of practical interest, and this paper develops a numerical 
method for calculating this admittance. The numerical method is based upon a 
Galerkin finite-element method and is described in the following section. 


NUMERICAL METHOD 

In the Galerkin finite-element method employed here, the transverse direction of 
the impedance tube (rectangular duct) is divided into a total of N-1 line elements, 
each with dimension b as shown in figure 2. These elements are considered 
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Figure 2.- Finite -element discretization and numbering system for grazing-flow 

impedance-tube test section. 



interconnected at joints which are called the nodes of the system. The nodes of the 
system are designated by the dark circles in the figure and are numbered appropri- 
ately. The number for each element is enclosed in a triangle and is identical with 
the lower node number of the element. In order to generalize and also simplify the 
formulation, it is convenient to define a local coordinate system for each element as 
shown in figure 3. Thus, global coordinates y^ and y I+1 for a typical element I 
of the global system are referred to as nodes 1 and 2, respectively, in this local 
coordinate system. Not only does this local coordinate system generalize and sim- 
plify the formulation but it will also facilitate the integration which is required 
to obtain the element equations. 


Global coordinate 
system 


Local coordinate 
system 



y 



y = y - y-,- 

b= v !+! " y l 


Figure 3.- Global and local coordinate systems for typical finite element. 


In this paper, a linear finite-element method is used. Thus, in each element 
the eigenfunction F(y) and Mach number M(y) are expressed in a form which ensures 
continuity of the functions F(y) and M(y) at the nodes joining two elements. In 
a typical element, these functions are expressed in the form 


F(y) = fj-fi^b) F I (I = 1, 2) (10) 


M(y) = f x (y,b) (I = 1 , 2) ( 11 ) 

in which repeated subscripts are to be summed over, Fj and M-j- are the values of 
the functions F(y) and M(y), respectively, at node I of the element, and 

f., (y*b) = 1 - (y/b) (12) 


f 2 (y,b) = y/b (13) 
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In following the method of Galerkin, the residual error E(y) in each element 
is multiplied by the basis function fj(y,b) and integrated over the element to 
obtain 


r 




/ E(y ) f (y) dy (1=1) 


/ E(y) f ;[ (y,b) dy = / 




(14) 


/ E(y ) f (y,b) dy 
0 * 


(I = 2) 


where 


2 1/ ~ \ 2 

E(y ) = ( 1 - K Ml F n + 2K M* F * + K 1 - K Mj 
x / x i V x 


~2 22,22 
K + m n /K L 
x 1 ' 


) 1 ( 


1 - KM F 


(15) 


Equation (14) expresses the desired averaging of the residual error, but it does 
not admit the influence of the boundary conditions* Integration by parts is now 
incorporated to reduce the order of the second-derivative term in equation (14) and 
to incorporate the boundary conditions. After integration by parts, equation (14) 
becomes 


/ E(y) f ] .(y) dy = / j^K^M'F’f i - ^1 - + K. 2 ^1 - 


- 1 


~ \ ~2 ( 2 2 , 2 2 \ 

- K Ml K + Im n /K L I 

' ' — 


Ff T dy 


+ F’(b) 


1 - K M(b) 

f (b,b) - F’(0) 

1 - K M(0) 

X 


X 


f x (0,b) 


(16) 


In regard to the first element, equation (8) is used to eliminate F’(0); therefore. 


F’ (0) 


1 - K M( 0) 
x 


f i (0 - b> - - iK K,x 


1 - K M( 0 ) 
x 


1 - K M( 0) 
x 


F(0) f ^ (0,b) 


(17) 
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Similarly, equation (9) is used to eliminate F' (b) from the last element so that 


F' (b) 


- - 


r ~ i 

n 

— — 

1 - K M(b) 

X 

£ l( b,b) = iKfl^ 

1 - K M(b) 

X 

_J 


1 - K M(b) 

X 


F(b) f x (b,b) 


(18) 


Substituting equations (10), (11), (17), and (18) into equation (14) leads to 
the element stiffness matrix. Thus, the contribution to the minimization of the 
residual error is expressed in matrix form: 


b 

/ E(y) f x (y,b) dy = [A q ]{$ q } 
0 


in which the element stiffness matrix [A q ] is 



[A q ] = 




the element vector {$ q } contains the values of the function F at nodes q and 
q+1 of the element 


«& q } 


Tr 



-q +1 


} 


and the superscript Tr denotes the vector transpose. 

Now, equation (19) gives the contribution to the minimization of the residual 
error due to element q. The contribution to the residual error due to elements 1 to 
N— 1 is obtained by requiring that at nodes where any two elements are connected, the 
value of the function F(y) at that node is the same for both elements. This 

requirement is achieved by adding the matrix elements a^ 1 and a q 2 in the global 

system for any two connected elements. The result is a set of linear matrix equa- 
tions, which are expressible in the form 


[A(p )]{«} = {0} 
w,u 


( 20 ) 
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where 


11 


1 2 


21 


22 


12 




21 


22 


1 2 


q+1 q 

a 22 3 11 + a 22 


(*> = (F,, F 2 , P 3 F^) 


N-1 

*21 


N-1 

22 


For a nontrivial solution to equation (20), the determinant (det) of the coeffi- 
cient matrix [A(8 )] must vanish: 

w,u 

detlAfp^^)] = {0} (21 ) 


Now, the power of the approach used here is that p w appears linearly in 

In ^ act ' °nly element in this coefficient matrix which depends on 

Pw,u is the element in the last row and column a^” 1 • The explicit dependence of 

this element on the unknown admittance can be determined. Thus, a N_1 is 
expressible in the form 22 


N-1 

22 


= b 


N,N 


- 8 d 
K w,u N,N 
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The coeffi- 


where N and d N N are independent of the unknown admittance P w u * 

cient matrix [A(B )] is decomposed as follows: 

w f u 

[A(p )] = [B] - p [D] (22) 

r w f u w,u 

N-1 

Thus, [B] is obtained from [A] by replacing a 22 by b N,N' whereas [D] has a single 
nonzero coefficient in its last row and column* 

Substituting equation (22) into (20) leads to the linear eigenvalue problem 


[B]{§> = P [D]{$} 
w,u 


(23) 


Because [B] is tridiagonal and [D] has a nonzero element in its last row and column 
only, equation (23) is fairly easy to solve. In this analysis, equation (23) is 
solved by using a special adaptation of the Gaussian elimination procedure that 
employs row pivoting and row scaling to maintain stability* This procedure trans- 
forms equation (23) into a new system of equations which are expressible in the form 


[i]{$} = 8 , ,[D]{®} 

W f u 


(24) 


b 1,1 b 1,2 b 1,3 


b 2,2 b 2,3 b 2,4 


[B] = 


b N-1,N-1 b N-1,N 


H,N 
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0 0 
0 0 


[D] = 


d N-1 ,N 
d N,N _ 


Since both [b] and [ 5 ] are upper triangular matrices, equation 
solved directly for the admittance p and eigenvector {$}. The 
this system reads WfU 


(23) can be 
last equation in 




(25) 


The solution to this equation is 


P 


w,u 



/d 


N,N 


(26) 


It is instructive to note that the value of p w given by equation (26) is 
unique and is the only value which makes the determinant of the coefficient matrix 
[A (p u ) ] vanish. (See eq. (21).) Further, the value of in equation (25) is 

set to unity and the remaining nodal values of the function F(y) are obtained from 
equation (24) by back substitution. 


RESULTS AND DISCUSSION 

In order to verify the numerical technique, results from the method are compared 
both with exact solutions that can be obtained for a constant mean-flow profile and 
with results of a Runge-Kutta integration technique for cases involving shear. To 
this end, only the lowest order cross mode is considered (m = 0), and continuity of 
particle velocity (n = 1 ) is employed as the boundary condition along the duct walls. 
Further, only 40 locations are employed in the finite-element discretization across 
the duct. The choice of this particular combination of parameters was felt to be 
sufficient to obtain confidence in the numerical technique. 
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Uniform Flow 


In the case of a constant mean-flow profile, the admittance (3 can be com- 
puted from the equation w,u 


P 


w,u 


\ 2 tan Kh + iK\P w ^1 - j -iK\^l - 


V 1 )" 


( ~ \ 

2n / „ \ 2n 

h - km) 

(l - KM) tan \h 

\ X / 

\ x / 


(27) 


where the eigenvalue \ is related to the decay rate by the equation 


X 2 = Hi -MK ' 2 

x, 


2 ( 2 2 2 2 \ 
K + I m n /K L I 


K 


Equation (27) is easily derived by finding the general solution to equation (7) 
for a constant mean-flow profile and then applying^the boundary conditions (eqs. (8) 
and (9)). Thus, for specified values of K, h, k x , L, m, and p w equa- 
tion (27) can be solved directly for the grazing-flow admittance. ' 

The values corresponding to a hard-wall duct are determined by setting both 

Pw,u an< ^ Pw,J£ to 0 in equation (27) and then solving this equation for K x . The 

K x values corresponding to the first- and third-order mode in a hard-wall duct have 
been computed from equation (27) with K and h both set to unity. These values 
are presented to two decimal places in table I for a range of constant mean-flow pro- 
files. The values of P w r K, h, M, and K x given in table I were used as 
input into equation (20) so that the admittance (3^ could be determined. In all 
cases, the admittance value determined from equation (20) agreed with the exact value 
of 0 to within four decimal places. 

The present method is also applicable to ducts for which the bottom wall is not 
hard (i.e., P w ^ ^ 0) but ^ as a finite value. Values of P w u computed from 
equation (27) for = 0 + I.Oi, with K and h both set* to unity, and for a 

range of M and K x 'values are presented in table II. The values of p w M, K x , 
K, and h given in table II were used as input into equation (20) and thus the 
admittance could be determined. In all cases, the admittance (3 determined from 

equation (20) agreed with those in table II to the number of decimal places shown in 
the table. 


Shear Flow 

When there is shear so that M’(y) is not 0, an exact expression for the admit- 
tance is apparently impossible. Therefore, the effort was made to evaluate the 
validity of the current numerical method in the presence of shear by comparing 
results to that of a Runge-Kutta integration technique (ref. 1). 
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TABLE I.- AXIAL PROPAGATION CONSTANT FOR FIRST- AND THIRD-ORDER 

MODE IN HARD-WALL DUCT 

[K = 1; 8 =8 = 0 + Oi: h = 11 

^w,u 


K for first 

X 

-order mode 

2.00 

+ 

O.OOi 

1 .67 

+ 

• 00i 

1 .43 

+ 

• OOi 

1.25 

+ 

• OOi 

1.11 

+ 

.OOi 

1.60 

+ 

.OOi 

.91 

+ 

.OOi 

.83 

+ 

.OOi 

.77 

+ 

.OOi 

.71 

+ 

.OOi 

.67 

+ 

.OOi 


K for third-order mode 
x 


0.67 + 7. 1 3i 
.48 + 6.75i 
.33 + 6.49i 
.21 + 6. 33i 
.10 + 6. 23i 
.00 + 6. 20i 
-.10 + 6. 23i 
-.21 + 6. 33i 
-.33 + 6. 49i 
-.48 + 6. 75i 
-.67 + 7. 1 3i 


TABLE II.- EXACT VALUE OF 8 „ FOR SOFT-WALL DUCT 

W j XX 

[K = 1; p = 0 + I.Oi; h = 11 

w,£ 


K 

X 

2.29 

+ 

1 • 83i 

1.99 

+ 

1.76i 

1.767 

+ 

1 .709i 

1.598 

+ 

1.675i 

1.465 


1 .655i 

1.356 

+ 

1.650i 

1.263 

+ 

1 .655i 

1.182 

+ 

1 • 675i 

1.107 

+ 

1 .709i 

1.035 

+ 

1 • 761 i 

.958 

+ 

1 . 835i 


Exact admittance, p 


0.519 - 
.662 - 
.860 - 
1.135 - 
1.512 - 
2.000 - 
2.545 - 
2.967 - 
3.041 - 
2.738 + 
2.255 + 


0. 967i 
1 . 1 51 i 
1 ,334i 
1 . 500i 
1 .614i 
1 .600i 
1 .349i 
.734i 
.061i 
.561i 
1 . 899i 













The k x values presented in table III were computed from a Runge-Kutta integra- 
tion technique identical to that used in reference 1 . Values of K x shown in the 


TABLE III.- AXIAL PROPAGATION CONSTANT FOR FIRST-ORDER MODE 
IN HARD-WALL DUCT FOR CONSTANT-GRADIENT SHEAR FLOW 

tK = 2%> P w,A = P w,u = 0 + 0i; h = 2 J 


M o 

** 

K for first-order mode 

X 

-0.3 

7.880 + O.OOOi 

-.1 

6.655 + . OOOi 

.1 

6.010 + .000i 

.3 

5.612 + .000i 


table correspond to the first-order mode in a hard-wall duct 

= ^w,u = 0 + 0*00i) with h = 2 and K = 2 it. A constant-gradient laminar 
shear flow* was employed for which 


M(y) 


"2M y/h 
o ' 


2M [1 - (y/h)] 
o 


(0 < y < h/2) 
(h/2 < y < h) 


(28) 


where M Q is the centerline Mach number. The values of h , K x , Mq, p w Zf 
and K shown in table III were input into equation (20) and thus the admittance 
Pw,u could be determined. In all cases, values of p^ determined from equa- 
tion (20) were 0 up to four decimal places. ' 


CONCLUDING REMARKS 

A new method has been developed for calculating acoustic-liner admittance in a 
rectangular duct with grazing flow. The approach is based upon a finite-element 
discretization of the acoustic field and a reposing of the unknown admittance value 
as a linear eigenvalue problem on the admittance value. Gaussian elimination is 
employed to solve this eigenvalue problem. Unlike existing methods that are based 
upon the integration of an ordinary differential equation and are thus limited to 
mean flows with one-dimensional boundary layers, the present method is extendable to 
mean flows with two-dimensional boundary layers as well. 

Admittance values determined from the method were compared with exact solutions 
obtained for a constant mean-flow profile, with excellent comparisons being obtained. 
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In the presence of shear, the results of the method were compared with results of a 
Runge-Kutta integration technique, with good comparisons again being obtained. The 
excellent results obtained from the method for the various flow profiles considered 
here motivate its extension to include grazing flows with two-dimensional boundary 
layers. Such an extension should now be formulated. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
May 22, 1984 
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